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Abstract 

Scanning  tunneling  microscopy  observations  of  a  reversible,  tip-induced  transition 
between  triangular-  and  star-shaped  networks  of  partial  dislocations  on  the  basal  <"0001) 
plane  of  highly  oriented  pyrolytic  graphite  (HOPG)  are  reported.  The  transition  between 
network  geometries  results  from  small  variations  (0.14  V)  in  the  tip-to-substrate  voltage 
bias  and  is  attributed  to  shear-induced  motion  of  partial  dislocations.  The  shear  stress 
required  for  the  transition  is  estimated  to  be  5  MPa.  Dislocation  motion  occurs  over 
distances  of  tens  of  nanometers  on  the  time  scale  of  several  minutes,  allowing  the  dynamics 
of  the  transition  between  the  two  networks  to  be  observed  in  real  time.  Atomically  resolved 
images  near  a  dislocation  show  distortions  of  the  atomic  lattice  that  are  consistent  with  glide 
of  the  surface  basal  plane.  Analysis  of  the  separation  distance  between  pairs  of  dislocations 
yields  a  surface  stacking  fault  energy  of  2. 1-4.5  mJ/m2. 
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I.  Introduction 


‘J 


\ 


Highly  oriented  pyrolytic  graphite  (HOPG)  has  been  used  extensively  as  a  substrate 
in  scanning  tunneling  mircoscopy  (STM)  investigations1-2'3  because  it  is  readily  cleaved  to 
yield  atomically  flat  stable  surfaces.  Although  Lipson  and  Stokes4  reported  in  1942  that 
natural  graphite  samples  are  comprised  of  5-15%  rhombohedral  graphite,  intermixed  with 
the  more  abundant  hexagonal  graphite,  STM  obseivations  of  the  rhombohedral  structure 
have  only  been  reported  once5.  The  percentage  of  rhombohedral  graphite  in  synthetic 
HOPG,  which  is  commonly  used  in  STM  investigations,  depends  on  the  annealing  and 
quenching  procedures  used  in  the  specimen  preparations,  and  is  typically  on  the  order  of  a 
few  percentile4'6-7. 

.Rhombohedral  graphite  differs  from  hexagonal  graphite  in  the  stacking  of  atomic 
p!  ar.es  parallel  to  the  c-axis.  Rhombohedral  graphite  has  an  abcabc  sucking  sequence 
while  hexagonal  graphite  has  an  abab  sucking  sequence  (Fig.  la).  Previous  transmission 
electron  microscopy  (TEM),  electron-  and  x-ray  diffraction  studies  have  demonstrated  that 
the  rhombohedral  phase  is  contained  within  sucking  faults  in  pyrolytic  graphite.  The 
sucking  faults  'are  generally  observed  in  TEM  as  networks  of  partial  dislocations  that 
separate  regions  of  faulted  (rhombohedral)  and  unfaulted  (hexagonal)  graphite.  By 
displacement  of  one  or  more  atomic  plane(s)  normal  to  the  c-axis,  the  transition  between  the 
rhombohedral  and  hexagonal  phases  within  the  network  may  be  represented  as 

ababa ...  t-tabcab...  (1). 

The  resulting  "glide"  of  an  atomic  plane,  shown  in  Fig.  1(b),  defines  the  partial 
dislocations  in  graphite  and  has  been  detailed  by  Delavignette  and  Amelinckx8, 
Williamson9,  and  Lipson  and  Stokes4.  An  in-plane  transition  between  the  two  phases, 
(which  are  bounded  by  the  partial  dislocation,  requires  distortion  of  the  atomic  arrangement 
and/or  cleavage  of  carbon-carbon  bonds.  However,  breakage  of  in-plane  carbon-carbon 
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bonds  is  unlikely  based  on  energetic  considerations.  Low  energy,  in-plane  distortions  of 
carbon-carbon  bonds  are  likely  to  occur  along  the  dislocation  for  small  distances  of  in-plane 
glide  (~  one  atomic  row),  although  they  are  probably  less  severe  than  suggested  in  Fig. 
1(b).  The  in-plane  atomic  bond  distances  (1.42  A)  and  interlayer  separation  (3.35  A)  of  the 
pure  hexagonal  and  rhombohedral  phases  removed  from  the  dislocation  are  identical. 

In  this  report,  STM  images  of  large-scale  triangular-  and  star-shaped  dislocation 
networks  that  separate  regions  of  hexagonal  and  rhombohedral  graphite  on  the  basal 
surface  of  an  HOPG  sample  are  reported.  A  reversible  transition  between  the  triangular- 
and  star-shaped  networks  is  observed  that  is  a  function  of  the  bias  voltage  between  the 
STM  tip  and  HOPG  substrate.  The  transition  results  in  the  interconversion  of  hexagonal 
and  rhombohedral  graphite,  indicating  the  reversibility  of  eq.  (1).  Atomically  resolved 
images  of  the  transition  region  yield  direct  structural  details  of  the  atomic  arrangement  along 
the  partial  dislocation.  The  energy  between  the  surface  dislocations  in  the  network  is  also 
determined  and  compared  to  bulk  values  obtained  from  previous  TEM  studies. 

II.  Experimental 

Tunneling  experiments  were  performed  using  a  commercial  STM  (Nanoscope  II, 
Digital  Instruments  Inc.,  Goleta,  CA).  Images  were  recorded  in  the  constant  current  mode 
using  a  mechanically  cut  Pt-Rh  (80-20%)  tip.  Tunneling  currents,  bias  voltages  (reported 
as  sample  versus  tip),  and  scan  rates  are  reported  in  each  figure  caption.  All  images,  except 
for  those  in  Fig.  3,  have  been  low  pass  filtered  once.  Images  in  Fig.  3  have  one  fast 
fourier  transform  applied  to  them.  Experiments  were  performed  in  air  cn  a  freshly  cleaved 
sample  of  HOPG  (Union  Carbide  grade,  B). 

III.  Results  and  Discussion 

Our  results  are  presented  in  the  following  order.  In  Section  IHA,  we  describe  the 
general  features  of  the  dislocation  network  investigated  in  this  study.  In  Section  IIIB, 
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atomically  resolved  images  of  the  faulted  structure  are  presented  which  provide  details  rf 
the  distortion  of  carbon-carbon  bonds  along  a  dislocation.  The  microscopically-  and 
atomically-resolved  images  presented  in  these  sections  are  consistent  with  the  structure 
shown  in  Fig.  1  and  with  the  notion  that  the  faulted  region  results  from  in-plane  glide  of  the 
basal  plane.  In  Section  IIIC,  images  of  reversible  STM-tip  induced  transitions  between 
triangular-  and  star-shaped  networks  are  presented  and  analyzed  in  terms  of  the  extension 
and  contraction  of  dislocation  nodes.  A  preliminary  analysis  of  the  tip-induced  transition 
dynamics  is  also  presented,  including  calculations  of  an  average  dislocation  velocity  from  a 
time-resolved  series  of  images.  In  Section  IIID,  we  calculate  the  energy  of  the  dislocation 
array  from  STM  images  and  esdmate  the  tip-induced  shear  forces  required  to  cause  the 
transition  between  network  geometries.  We  show  that  the  required  shear  forces  are  of  the 
same  order  of  magnitude  as  tip-substrate  forces  required  to  image  HOPG  samples. 

A.  Dislocation  Network. 

Fig.  2(a)  shows  an  STM  image  of  a  triangular-shaped  network  of  partial 
dislocations  on  the  HOPG  surface.  On  average,  the  bases  and  heights  of  the  triangles  were 
75  and  1 10  nm,  respectively.  Triangular  arrays  similar  to  that  shown  in  Fig.  2  have  been 
previously  observed  in  STM5  and  TEM8-9  investigations.  The  triangular  network 
encompassed  a  total  area  estimated  to  be  in  excess  of  4  pm2. 

One  edge  of  the  network  was  examined  in  detail.  Along  this  edge,  a  row  of 
elliptical  dislocation  loops  defined  the  boundary  between  one  comer  of  each  triangle  and  a 
long  straight  ribbon.  Fig.  2(b).  The  dislocation  loops  are  -50  nm  wide  and  -47  nm  long, 
and  have  an  apparent  displacement  of  3.9  A  above  the  surrounding  surface.  The  central 
depression  defined  by  the  loop  is  -2.9  A  deep.  Elliptical  and  quasi-hexagonal  dislocation 
loops  have  been  previously  identified  on  the  graphite  basal  plane  by  Delavignette  and 
Amelinckx8  and  are  believed  to  pin  dislocations  associated  with  in-plane  glide,  consistent 
with  the  termination  of  the  triangular  network.  The  ribbons  emerging  from  the  triangular 
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array  are  comprised  of  two  partial  dislocations  (vide  infra)  separated  by  12-25  nm;  the 
center-to-center  distance  between  ribbons  is  equal  to  the  base  width  of  a  triangle  within  the 
network  (75  nm).  The  partial  dislocations  form  a  large  network  of  parallel  ribbons  (at  least 
2.5  x  2.5  }im,  Fig.  2(c))  and  appear  to  terminate  at  a  poorly  resolved  defect  (not  shown) 
that  is  probably  associated  with  plastic  deformation  during  cleavage  of  the  sample.  This 
conclusion  is  supported  by  previous  TEM  observations  that  suggest  that  the  number 
density  of  stacking  fault  arrays  in  pyrolytic  graphite  significantly  increases  with  mechanical 
crushing  of  the  sample  and  that  dislocations  end  at  the  edge  of  a  flake,  cleavage  step  or  low 
angle  boundary9. 

Williamson9  attributed  the  diffraction  contrast  in  TEM  images  of  triangular-shaped 
arrays,  which  are  nearly  identical  to  the  array  in  our  STM  images  (Fig.  2),  to  rhombohedral 
stacking  faults  enclosed  between  partial  dislocations  in  the  hexagonal  structure.  The  glide 
of  atomic  planes  observed  in  HOPG  by  TEM  is  generally  in  the  form  of  arrays  of  partial 
dislocations  which  can  combine  into  a  perfect  basal  screw  dislocation  following  the  general 
reaction 


l/3af  1100]  +  l/3a[1010]  =  l/3a[2110]  (2). 

The  partials  describe  the  translocation  needed  to  convert  a  'b'  layer  to  a  'c'  layer,  eq.  (1), 
and  form  the  boundary  between  the  hexagonal  (unfaulted)  and  the  rhombohedral  (faulted) 
phases.  Based  on  the  structural  assignments  from  previous  TEM  and  electron  diffraction 
analyses8*9,  and  the  geometrical  sl-nilarity  of  the  network  array  in  both  STM  and  TEM 
images,  we  have  previously  assigned  the  lighter  triangles  in  the  STM  images  to  the 
rhombohedral  crystal  structure  and  the  darker  triangles  to  the  hexagonal  crystal  structure5. 
Similarily,  the  rhombohedral  and  hexagonal  structures  are  assigned  to  the  bright  and  dark 
regions,  respectively,  of  the  ribbons  extending  from  the  triangular  network  edge. 


The  observed  contrast  in  the  image  of  the  triangular-shaped  network  is  dependent 
on  the  STM  sample  bias  voltage  and  tunneling  current.  In  the  image  shown  in  Fig.  2(a), 
which  was  recorded  at  a  bias  voltage,  Vt  (substrate  vs.  tip),  of  0.24  V  and  set-point 
current,  It,  of  0.37  nA,  the  lighter  contrast  triangles  correspond  to  an  apparent  displacement 
of  1.7  A  above  the  darker  triangles.  The  intersection  of  three  triangles  is  defined  by  a 
brighter  circular  region  (diameter  -  20  nm)  which  is  displaced  4.3  A  above  the  dark 
triangular  regions.  When  the  STM  settings  were  changed  to  Vt  =  0.53  V  and  It  =  2.4  nA, 
the  contrast  between  the  light  and  dark  triangles  disappears  and  the  intersection  points 
appeared  displaced  by  1.2  A  above  the  surface.  At  Vt  =  -0.22  V  and  It  =  0.37  nA,  the 
contrast  difference  between  the  dark  and  light  triangles  reduced  to  0.8  A  and  the 
intersection  points  disappeared.  We  note  that  the  interplanar  spacings  of  hexagonal  and 
rhombohedral  graphite  are  identical,  and  that  the  triangular  network  structure  is  consistent 
with  the  assumption  that  glide  is  contained  within  the  basal  plane.  Thus,  the  observed 
contrast  cannot  result  from  an  actual  topographical  z-displacement  between  the  faulted  and 
unfaulted  regions.  A  similar  contrast  depended*  on  Vt  and  It  was  also  noted  in  previous 
work5  and  discussed  in  terms  of  the  differences  in  electronic  and  elastic  properties  of 
hexagonal  and  rhombohedral  graphite.  However,  no  unambiguous  explanation  of  the 
observed  contrast  has  yet  been  forwarded. 

A  recent  publication  by  Garbarz  et  al.10  has  associated  a  superlattice  hexagonal 
network  (-35  A  periodicity)  with  an  array  of  basal  dislocations.  Similar  large-scale 
hexagonal  structures  with  various  periodicities  have  been  observed  in  our  laboratory  and  by 
others1 1  and  identified  as  a  local  misorientation  of  the  top  basal  plane  giving  rise  to  a  simple 
moirti  pattern.  However,  we  have  never  observed  these  structures  near  a  large-scale 
dislocation  network,  such  as  that  shown  in  Fig.  2. 

The  STM  observations  that  we  describe  in  detail  below  are  based  soley  on 
continuous  studies  of  the  single  dislocation  network  shown  in  Fig.  2  over  a  period  of  two 
days.  While  we  have  previously  reported  STM  images  of  similar  structures5,  well-defined 


network  arrays  of  this  magnitude  are  infrequently  encountered  during  STM  imaging. 
Experimental  n  .easurements  on  the  dislocation  terminated  with  a  tip  crash  which  resulted  in 
a  marked  decrease  in  image  resolution.  We  were  unable  to  relocate  the  network  array  using 
new  atomically  sharp  tips. 

B.  Atomically-Resolved  Structure  of  a  Partial  Dislocation. 

Fig.  3(a)  shows  an  atomically  resolved  image  near  the  center  of  a  dark  (unfaulted), 
triangular  region.  The  atomic  arrangement  and  spacing  is  consistent  with  the  reported 
hexagonal  structure  of  graphite,  i.e.,  the  periodicities12  are  2.6  A  and  the  corresponding 
fast  Fourier  transform  (FFT)  of  the  image  (inset)  shows  six  equally  spaced  points  with 
hexagonal  symmetry.  However,  the  atomic  corrguation  is  0.8  A,  larger  than  expected. 
Large  corrugation  amplitudes  have  previously  been  observed  on  hexagonal  graphite  and 
have  been  attributed  to  elastic  interactions  of  the  tip  with  the  sample  surface13*14.  In 
regions  'ar  removed  from  the  dislocation  network,  atomically  resolved  images  cf  the 
graphite  surface  show  a  hexagonal  symmetry  with  periodicities  of  2.6  A  and  a  corrugation 
amplitude  of  0.3  A. 

The  observed  periodicity  and  symmetry  in  the  STM  images  is  known  to  result  from 
an  inequivalence  of  the  electron  charge  density  between  surface  atoms  on  the  basal  plane  of 
hexagonal  graphite.  The  tunneling  current  density  above  surface  atoms  with  neighboring 
atoms  located  directly  below  in  the  second  layer  (and  all  subsequent  layers),  i.e.,  a-sites,  is 
significantly  lower  than  for  surface  atoms  which  do  not  have  neighboring  atoms  in  the 
second  layer,  i.e.,  P-sites  (Fig.  1(a)).  Thus,  the  bright  spots  observed  in  the  STM  images 
of  hexagonal  graphite  correspond  to  regions  of  high  current  density  above  P-site  atoms. 

In  bulk  rhombohedral  graphite,  the  inequivalence  between  in-plane  neighboring 
atoms  vanishes  since  all  atoms  in  the  bulk  have  one  nearest  neighbor  atom  located  directly 
in  an  adjacent  plane.  However,  the  calculations  of  Toman6k  and  Louie15  suggest  that  the 
asymmetry  in  the  charge  density  between  neighboring  in-plane  atoms  persists  at  the  surface 
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of  rhombohedral  graphite  due  to  a  surface  state  near  the  Ft.  level  that  is  localized  on 
surface  atoms  which  do  not  have  neighboring  atoms  in  the  second  layer  (Fig.  1(a)).  Thus, 
a  similar  corrguauon  and  six-fold  symmetry,  as  shown  in  Fig.  3(a)  for  hexagonal  graphite, 
would  also  be  expected  for  rhombohedral  graphite.  Unfortunately,  we  were  not  able  to 
ootain  a  clear,  atomically-resolved  image  of  a  region  within  one  of  the  bright  triangles  to 
test  this  prediction.  However,  Fig.  3(b)  shows  an  atomically -resolved  image  of  a 
boundary  between  a  light  and  dark  triangle.  (The  image  has  oeen  filtered  using  a  FFT  to 
improve  contrast  of  the  atomic  scale  corrugation;  this  filtering  reduces  the  contrast  defining 
the  boundary  between  the  bright  and  dark  triangles  of  the  network.  The  periodic  atomic- 
scale  strucutre  observed  in  the  image  prior  to  FFT  filtering  was  not  significantly  affected  by 
this  process)  The  atomic  arrangement  in  the  upper  region  of  the  image  (corresponding  to 
the  faulted  rhombohedral  structure)  shows  a  distorted  hexagonal  arrangement  with 
periodicities  of  2.4  A  and  2.8  A.  The  atomic  corrugation  is  0.5  A  near  the  dislocation  line 
drawn  in  the  figure  and  decays  to  0.3  A  away  from  the  line  in  either  direction.  The  FFT  of 
the  image  (inset)  shows  four  points  in  a  rectangular  arrangement,  consistent  with 
rectangular  electron  diffraction  patterns  previously  observed  near  dislocations  on  pyrolytic 
graphite16.  A  distortion  in  the  hexagonal  lattice  can  be  seen  by  following  the  line  drawn 
from  the  top  left  to  the  bottom  right  of  the  image  at  an  angle  of  62°  to  the  direction  of  the 
dislocation  (which  is  in  the  <1100>  direction).  The  atomic  rows  along  the  <1010> 
direction  shift  a  distance  equivalent  to  one  unit  cell  (~2.7  A  measured;  2  46  A  theory)  in 
crossing  the  dislocation,  consistent  with  the  model  structure  shown  in  Fig.  1(b).  The  shift 
occurs  over  a  distance  of  12  hexagonal  unit  cells,  which  we  interpret  as  defining  the  width 
cf  the  dislocation.  This  finding  is  in  agreement  with  previous  estimates  of  the  widths  of 
partial  dislocations  in  graphite  (-20-30  unit  cells17).  Thermal  and  electronic  drift  in  these 
measurements  was  negligible  (see,  for  example,  Fig.3(a)),  ruling  out  the  possibility  that  the 
distortion  of  the  hexagonal  structure  is  a  measurement  artifact. 
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Cose  inspection  of  Fig.  3(b)  shows  that  the  hexagonal  structure,  albeit  distorted,  is 
maintained  throughout  the  dislocation  region.  This  finding,  along  with  the  observation  that 
the  contrast  remains  approximately  uniform  across  this  boundary,  suggests  that  the  charge 
density  asymmetry  between  neighboring  surface  atoms  occurs  within  the  dislocation 
region,  as  well  as  on  the  surface  of  the  hexagonal  (observed)  and  rhombohedral  (predicted) 
phases.  However,  details  of  the  atomic  structure  along  the  dislocation  are  not  sufficiently 
well  understood  to  analyze  this  finding  in  terms  of  electronic  arguments. 


C.  Tip-Induced  Dislocation  Motion  and  Network  Transitions. 

I 

A  reversible  transformation  of  the  triangular-shaped  network  to  a  star-shaped 
network  was  observed  while  varying  the  bias  voltage  between  0.24  and  0.1  V,  with  the 
tunneling  current  maintained  at  2.4  nA.  Fig.  4(a)  shows  a  500  nm  x  500  nm  stable  image 
of  the  dislocation  network,  shown  previously  in  Fig.  1(a),  several  minutes  after  the  bias 
had  been  decreased  from  0.24  V  to  0  1  V.  The  network  now  appears  to  consist  of 

alternating  small  ard  large  triangles  oriented  at  180°  from  one  another  in  the  form  of  a 

| 

distorted  star.  The  sum  of  the  lengths  of  the  bases  of  small  and  large  triangles  (31  and  45 

I 

nm,  respectively)  is  equivalent  to  the  length  of  the  base  of  an  individual  triangle  (75  nm)  in 
the  network  in  Fig.  2(a).  A  similar  equivalence  exists  between  the  heights  of  the  triangles. 

The  sequence  of  images  in  Fig.  4(b-d)  demonstrate  the  reversibility  of  the  transition 
between  the  triangular-  and  star-shaped  networks.  Fig.  4(b)  is  a  enlargement  of  the  image 
in  Fig.  4(a)  showing  the  star-shaped  network  at  a  bias  voltage  of  0.1  V.  The  image  in  Fig. 
4(c)  was  recorded  at  -1  minute  following  a  change  in  bias  from  0.1  V  to  0.24  V.  Each 
triangle  is  observed  to  have  a  lighter  contrast  trapezoidal-shaped  band  -14  nm  in  width 
along  its  base  (enlarged  in  Fig.  5).  In  fact,  during  the  course  of  the  transition  from  the  star- 
to  triangular-shaped  network,  the  width  of  the  band  continuously  decreased  from  a 
maximum  value  of  -30  nm.  Bands  from  adjacent  triangles  overlap  and  the  bright 
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intersection  pc  ave  not  yet  reappeared.  The  image  in  Fig.  4(d)  was  taken  ~1  minute 
after  the  image  in  Fig.  <4vc)  with  the  bias  stili  at  0.24  V.  The  band  at  the  base  of  the 
triangles  has  decreased  to  9  nm,  a  few  of  the  mangles  have  joined  vertices  forming  the 
bright  intersection  points,  and  the  heights  and  bases  of  the  triangles  have  almost  attained 
their  original  values.  Continued  scanning  at  a  bias  voltage  of  0.24  V  for  several  more 
minutes  caused  a  complete  return  of  the  original  triangular-shaped  network  shown  in 
Fig.2(a).  This  reversible  transition  between  the  triangular-  and  star-shaped  network  was 
reproducible  during  continuous  imaging  over  a  period  of  two  days  with  no  degradation  of 
the  network  structures 

The  mechanism  for  the  transition  between  network  geometries  undoubtedly 
involves  the  movement  of  dislocations  within  the  network.  Fig.  6(a)  is  a  schematic  of  the 
triangular-shaped  network  showing  an  array  of  alternating  shaded  and  unshaded  triangles. 
The  unshaded  triangles  (extended  P  nodes18)  denote  the  rhombohedral  (faulted)  phase  of 
graphite  while  the  shaded  triangles  denote  regions  of  hexagonal  graphite.  The  boundaries 
enclosing  the  triangles  are  partial  dislocations.  At  the  intersection  of  the  dislocations, 
corresponding  to  the  bright  spots  in  Fig.  2(a),  the  dislocations  are  contracted  to  form  a  K 
node.  In  the  star-shaped  network.  Fig.  6(c),  both  the  K  and  ?  nodes  are  fully  extended 
and  enclose  rhombohedral  graphite. 

The  transition  from  the  triangular-shaped  network  to  the  star-shaped  network 
requires  expansion  of  the  contracted  dislocations  within  the  K  nodes,  as  indicated  in  Fig.  6. 
The  transition  can  be  viewed  as  resulting  from  the  concerted  motion  of  parallel  dislocations 
along  one  of  the  three  equivalent  <1010>  directions.  Indeed,  we  were  able  to  verify  from 
STM  images  that  motion  of  only  one  set  of  parallel  partial  dislocations  occured  during  the 
transition.  Not  coincidentally,  motion  occured  only  for  dislocations  oriented  nearly  at  right 
angles  to  the  direction  in  which  the  tip  was  scanning  (consistent  with  our  analysis  in 
Section  HID  that  motion  results  from  tip-induced  shear  stress).  We  note  here,  however, 
that  although  the  shear  forces  due  to  scanning  only  act  on  a  small  fraction  of  the  dislocation 
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network  at  any  one  instant  during  imaging,  the  transition  between  the  network  geometries 
results  from  the  concerted  motion  of  numerous  partials  (for  example,  the  concerted 
translation  of  3  parallel  partials  is  observed  in  the  network  transidon  shown  in  Fig.4(b-d)). 
The  translation  of  a  set  of  parallel  dislocations  along  only  one  of  the  three  equivalent 
<1010>  directions  is  shown  schematically  in  Fig.  6(b). 

The  stability  of  the  fully  developed  triangular-  and  star-shaped  networks  during 
prolonged  STM  scanning  suggests  that  these  structures  represent  metastable  states.  We 
believe  that  the  structure  shown  in  Fig.  6(b)  depicts  an  unstable  transition  state  captured  in 
the  image  shown  in  Fig.  4(c)  and  enlarged  in  Fig.  5.  The  schematic  structure  of  the 
intermediate  shows  that  motion  of  only  one  parallel  set  of  dislocations  is  consistent  with  the 
observation  of  a  quasi-trape~oidal  region  at  the  base  of  the  triangle.  As  the  partials  move 
across  the  surface,  the  width  of  this  trapezoidal  region  is  expected  to  change,  in  agreement 
with  the  time-resolved  images  of  the  network  transition. 

The  extent  to  which  the  triangular-to-star  transition  occurs  outside  of  the  imaged 
area  (following  switching  of  the  bias  between  0.24  and  0.1  V)  was  not  investigated. 
However,  the  ability  to  completely  and  rapidly  cause  the  reverse  transition  suggests  that  the 
transition  may  be  contained  within  a  small  local  region.  The  reversibility  of  the  transition 
also  suggests  that  one  of  the  two  metastable  states  (triangles  vs.  stars)  results  from  a 
cancellation  of  shear  forces  and  restoring  forces  that  result  from  a  local  distortion  of  the 
network.  We  speculate  that  pinning  of  the  dislocations  at  the  network  edgc(s)  by 
dislocation  loops,  such  as  those  observed  in  Fig.  2(b)  (as  well  as  other  surface  features), 
might  prevent  motion  of  the  dislocations  near  the  network  boundary.  This  would  require 
the  direction  of  the  dislocation  to  deviate  from  the  low  energy  <1010>  direction  near  the 
boundary,  raising  the  energy  of  the  network.  When  the  shear  force  is  reduced,  as  a  result 
of  changing  the  bias  voltage  (e.g.,  from  0.1  to  0.24  V),  the  energy  stored  in  the  network 
by  distortion  of  the  partials,  would  necessarily  be  sufficient  to  drive  the  system  back  into 
the  more  stable  state. 


We  were  able  to  estimate  average  velccitcs  of  dislocation  motion  in  the  above 
experiments.  From  Fig.  4(b-d),  for  example,  the  dislocation  velocity  can  be  approximated 
by  the  ratio  of  the  distance  the  triangles  elongate  to  fus“  at  their  vertices,  relative  to  the  time 
required  for  the  complete  transition  between  the  triangular-  and  star-shaped  networks  (-2 
minutes).  From  our  observations,  we  estimate  an  average  value  of  10'1®  m/sec,  which  is 
considerably  smaller  than  values  reported  in  the  literature  for  a  variety  of  materials  (Ifr8  to 
103  m/sec)19.  However,  in  macroscopic  measurements  of  dislocation  velocity,  a  shear 
stress  is  applied  across  the  entire  sample  such  that  the  resulting  shear  forces  act 
(presumably  uniformly)  along  the  entire  length  of  the  dislocation.  Obviously,  in  STM 
experiments,  at  any  one  instant,  the  tip-induced  shear  stress  is  highly  localized  over  a  short 
length  of  the  dislocation.  Thus,  it  is  not  unexpected  that  the  time  averaged  dislocation 
motion  resulting  from  a  series  of  short-time  pulses  over  the  dislocation  is  correspondingly 
smaller.  We  did  not  investigate  the  distinct  possibility  that  the  dislocation  velocity  may 
depend  on  the  scan  rate  as  well  as  the  applied  bias. 

D.  Stacking  Fault  Energy  and  Shear  Stress  for  Dislocation  Motion. 

Shear  stresses  resulting  from  interactions  between  the  tip  and  graphite  surface  are 
responsible  for  dislocation  motion.  Graphite  shear  deformation  ,  without  cleavage  of  the 
basal  plane,  can  only  occur  if  there  is  a  large  compressive  stress  normal  to  the  basal 
plane20.  In  addition,  the  rhomoohedral  phase  is  only  produced  when  the  hexagonal 
structure  is  under  severe  shear.  Thus,  the  triangle-to-star  shaped  network  transition 
requires  that  the  tip-sample  interaction  is  associated  with  both  a  normal  compressive  and 
shear  stress.  Evidence  to  support  the  presence  of  tip  induced  forces  is  evident  in  Fig.  2b 
and  4a,  where  the  dislocations  at  the  bases  of  the  triangles  are  bowed  out  normal  to  the 
direction  of  scanning. 

If  dislocation  motion  is  due  to  shear  stresses  resulting  from  tip  scanning,  then 
changing  the  tip-to-sample  distance  should  reveal  this  phenomenon.  Indeed,  the  voltage 
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dependent  transition  between  the  triangular-  and  star-shaped  networks  can  be  interpreted  in 
terms  of  the  dependence  of  the  tip-to-sample  separation,  s,  on  the  voltage.  Using  the 
tunneling  equation21 


I  o  Vexp(-2tCoS) 


(3) 


where  I  is  the  tunneling  current,  s  is  the  tip-to-sample  separation,  V  is  the  sample  bias 
voltage,  and  2k<j  (A)  =  1.025V$eff  (4»eff(eV)  is  the  effective  barrier  height  expressed  in 
terms  of  the  work  functions  of  the  two  electrode  materials  i.e.,  +  <(>hopg)/2)- 

From  eq.  (3),  the  change  in  tip-to-sample  separation,  As,  at  two  different  voltages  can  be 
estimated  if  tq,  is  known.  Since  the  value  of  the  barrier  height  varies  with  the  tip-to-sample 
separation  and  may  also  be  affected  by  the  presence  of  contaminants,  it  is  not  possible  to 
determine  an  absolute  value  for  however,  if  d>eff  is  assumed  to  remain  constant  in  the 
bias  voltage  ranges  used  (0.24  V  to  0.10  V),  then  the  absolute  change  in  s  can  be 
calculated.  For  example,  assuming  <peff  -  4  eV  and  a  constant  tunneling  current  of  2.4  nA, 
reducing  the  voltage  from  0.24  V  to  0.1  V  corresponds  to  a  -0.5  A  decrease  in  the  tip-to- 
sample  separation.  A  larger  value  in  As  would  be  expected  if  <{>cff  were  reduced  by  the 
presence  of  a  contamination  layer22.  We  speculate  that  this  bias-dependent  variation  in  tip- 

to-substrate  distance  influences  the  magnitude  of  the  shear  stresses  acting  on  the  surface 

1 

while  imaging. 

The  shear  force  necessary  for  dislocation  movement  can  be  calculated  from  the 
stacking  fault  energy,  Ysfe-  Numerical  values  of  Ysfe  can  be  obtained  from  the  separation 
distance  between  pairs  of  parallel  partial  dislocations,  such  as  those  shown  in  the  ribbon 
network  in  Fig.  2c.  Figure  7  shows  an  enlarged  view  of  three  sets  of  parallel  ribbons, 
each  ribbon  being  formed  from  a  pair  of  partial  dislocations  whose  separation  ranges  from 
12  to  25  nm.  The  stacking  fault  energy,  Ysfe.  can  be  calculated  from  eq.  (4)23 


YSFE  =  [(2-3Veff)/(l-Veff)]Meflb2/8Jld  (4) 

where  d  is  the  observed  spacing  between  the  partial  dislocations,  and  b  is  the  modulus  of 
the  Shockley  partial  Burgers  vector  (0.142  nm).  In  conjunction  with  anisotropic  elasticity 
theory  (Appendix  I),  the  elastic  constants  for  graphite  are:  Poissons  ratio,  vcff  =  0.34  and 
the  shear  modulus,  Peff  =  45.5  x  109  N/m2.  Using  these  values  and  the  range  of  measured 
separation  distances  (d  =  12-25  nm),  Ysfe  is  estimated  to  be  between  2.1  and  4.5  mJ/m2. 
Analysis  of  previously  published  STM  images5  of  a  similar  ribbon  network  yields  Ysfe  = 
1.5  -  4.7  mJ/m2 

The  average  value  of  Ysfe  obtained  from  our  STM  images  is  slightly  larger  than 
values  obtained  from  electron  microscopy  images24  (1.5  mJ/m2,  see  Ref.  25).  In  contrast 
to  electron  microscopy  experiments  where  the  observed  partial  dislocations  may  be  tens  to 
hundreds  of  angstroms  below  the  surface,  only  dislocations  at  or  within  the  first  few  layers 
of  the  surface  are  imaged  in  STM.  According  to  Read26  and  others27-28  the  width  of  the 
dislocation  ribbons  should  decrease  as  they  approach  the  surface  due  to  the  reduced  elastic 
replusion  between  the  dislocations.  Thus,  the  energies  of  the  surface  stacking  faults 
calculated  from  STM  observations  may  be  appreciably  larger  than  the  bulk  values  obtained 
by  electron  microscopy.  However,  it  has  been  theorized  that  the  stacking  fault  energy 
should  not  change  appreciably  for  very  anisotropic  materials,  such  as  graphite,  unless  the 
fault  is  located  very  close  to  the  surface30.  Correcting  for  the  surface  proximity  effect 
(Appendix  II),  we  obtained  Ysfe  *  0.14  -  2.6  mJ/m2  from  our  STM  images,  in  good 
agreement  with  electron  microscopy  results. 

To  estimate  the  force  required  for  dislocation  motion  we  consider  the  forces  acting 
on  the  partial  dislocations  enclosing  the  stacking  fault  ribbons.  An  elastic  repulsion  force 

between  the  partial  dislocations,  Fel.  tends  to  separate  the  partial  dislocations  while  an 
opposing  stacking  fault  force,  FyS(.,  tends  to  draw  the  partial  dislocations  together.  The 

stacking  fault  energy,  Ysfe>  has  units  of  energy  per  unit  area  (Nm/m2)  but  can  also  be 
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thought  of  as  a  force  per  unit  length  of  a  dislocation  (N/m)  which  draws  the  partial 
dislocations  together.  The  elastic  repulsion  force  also  equals  the  shear  force  per  unit 
length,  xb,  multiplied  by  the  length  of  the  dislocation  which  is  affected  by  the  force,  Leff. 
From  Kelly  and  Groves29,  the  stacking  fault  force  is  also  equal  to  Ysfe^  giving 

FYsf  (N)  =  FEL  (N)  =  TbLeff  (N)  =  ysFE  d  (N).  (5) 

The  field  around  the  STM  tip  can  be  thought  of  as  a  repulsive  force  which  destabilizes  the 
partial  dislocations  from  their  equilibrium  separation.  Thus,  the  tip  acting  on  a  length  of  a 
dislocation,  Leff,  causing  it  to  move,  would  exert  an  average  shear  stress  of 

X  =  YSFE  tVbLeff  (6) 

As  a  zero-oraer  estimate  of  T,  we  assume  that  leff  is  equivalent  to  the  length  of  the  base  of 
the  triangle  between  two  bright  points  (75  nm),  yielding  x  =  5  MPa.  (For  comparison, 
macroscopic  indentation  experiments  indicate  that  forces  on  the  order  of  ~1  GPa  damage 
the  surface  of  graphite30).  Mate  et  al.  have  shown  that  forces  between  the  tip  and  substrate 
on  the  order  of  10^  to  10-7  N  (-50  to  100  MPa  for  their  assumed  80  nm  tip  contact  radius) 
were  required  to  obtain  STM  images31.  In  addition,  Gilman32  has  estimated  that 
dislocation  motion  can  be  activated  at  room  tempera  tun*  it x  >  p/100  or  500  MPa.  From 
this  analysis,  the  shear  force  provided  by  the  STM  tip  is  both  less  than  the  force  Mate  et  al. 
found  necessary  for  tunneling  and  less  than  the  force  required  for  dislocation  motion.  The 
calculated  value  for  x  using  Eq.  (6)  can  also  be  compared  to  the  macroscopic  average 
critically  resolved  shear  stress  for  the  basal  plane  of  annealed  HOPG  which  is  -0.6  MPa33, 
or  -lOx  less  than  the  forces  estimated  from  our  STM  experiment.  At  best,  the  calculation 
in  Eq.  (6)  is  only  a  crude  estimate  as  the  tip  is  acting  on  a  small  portion  of  the  dislocation 
line  (which  may  be  less  than  75  nm)  at  one  time.  The  value  of  5  MPa  could  be  considered 
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a  lower  bound  for  x.  To  estimate  an  upper  bound  for  x,  one  can  take  Leff  =  2  nm  (which 
assumes  that  only  the  section  of  the  dislocation  directly  under  the  tip  is  moving),  and  obtain 
a  value  for  x  of  ~200  MPa,  comparable  to  or  larger  than  the  abovementioned  estimates  of 
the  force  necessary  for  dislocation  motion  and  STM  imaging.  We  conclude  that  the  shear 
forces  resulting  from  the  tip  are  sufficient  to  cause  the  observed  transitions  between 
triangle-  and  star-shaped  networks. 

IV.  Summary 

We  have  reported  STM  observations  of  a  reversible,  tip-induced  transition  between 

two  different  network  geometries  of  partial  dislocations  on  the  basal  plane  of  HOPG.  The 

shear  stress  required  for  dislocation  motion  is  comparable  to  previous  estimates  of  the  shear 

stress  in  STM  experiments.  Atomically  resolved  images  near  a  dislocation  show 

distortions  in  the  lattice  which  are  attributed  to  glide  of  the  dislocation  on  the  basal  plane. 

In  addition,  dislocation  motion  was  observed  on  the  time  scale  of  the  STM  experiment, 

allowing  an  estimate  of  the  dislocation  velocity.  The  capability  of  STM  to  image  and 

control  dislocation  motion  in  real  time  may  provide  a  new  approach  to  probe,  at  the 

nanometer  scale,  the  effects  of  various  treatments  (e.g.,  quenching  and  annealing)  on  the 

material’s  properties.  For  instance,  it  should  be  possible  to  observe  how  dislocation 

% 

pinning  by  impurities  or  the  interaction  of  dislocations  with  surface  adsorbates  affects 
dislocation  motion  and  thus  the  elastic  and  plastic  properties  of  the  material. 
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APPENDIX  I:  Calculation  of  ySFE 

The  width,  d  (=  12-25  nm  from  Fig.  7),  of  a  ribbon  lying  in  the  basal  plane  of  an 
anisotropic  hexagonal  crystal  is  given  by26, 

d  =  do[l-(2veff/(2-veff))cos2<p]  (I.l) 

with 

do  =  [(2-veff)/(  l-veff)]peffb2/8rtYSFE  (1.2) 

where  veff  is  Poisson's  rado,  9  is  the  angle  between  the  total  Burgers  vector  and  the 
direction  of  the  ribbon,  do  is  the  equilibrium  spacing  between  the  partials,  b  is  the  modulus 
of  the  Shockley  partial  Burgers  vector  (0.142  nm),  and  peff  is  the  shear  modulus.  For  a 
pure  screw  disloction,  9  is  zero  since  the  Burgers  vector  is  parallel  to  the  direction  of  the 
ribbon  and  so  d  =  do[l-(2veff/(2-veff)].  In  an  anisotropic  medium,  Peff  and  veff  must  be 
calculated  according  to26 


and 

where 

and 


Meff  =  [C44/2(Ci  1-Ci2>]1/2 

Heff /(1-Vcff)  =  (C11C33  -  Ci32)/C33(26i  +  §2)1^2 

5i  =  (C11/C33)1/2 


(1.3) 

(1.4) 

(1.5) 


§2  =  (CUC33  -  C132  -  2C13C44VC33C44 


(1.6) 


In  eq.’s  (I-3-I.6),  c&  are  the  elastic  constants  for  a  system  with  the  *3  axis  along  the 
hexagonal  axis.  Values  for  c,k,  which  are  taken  from  Kelly34  and  listed  below,  are  given 
as  Cik  =  Cfr/101 1  dyn  cm'2: 

ci  1  =  106  ci2=18  C13  =  1.5  C33  =  3.65  C44  =  0.45. 

The  Cauchy  relations  for  central  forces,  i.e.,  C13  =  C44  and  C12  =  l/2(cn  -  C12),  are  not 
obeyed.  The  value  for  C44  shown  above  is  for  a  perfect  crystal  with  the  glissile  basal 
dislocations  pinned  out.  The  effects  of  twist  boundaries,  cross-links  between  adjacent 
layer  planes,  and  specific  dislocation  configurations  on  the  value  of  C44  are  not  known. 

Using  the  above  equations  (1. 1-1.6)  in  conjunction  with  the  values  for  Cfc,  jieff  = 
45.5  x  10-9  N/m2,  vcff  =  0.34,  and  YSFE  =  2.1  -  4.5  mJ/m2 
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APPENDIX  II:  Correction  of  YSFE  for  Surface  Proximity  Effects 

Calculation  of  Ysfe  IS  based  on  the  separation  distance,  d,  between  the  partial 
dislocations  in  a  ribbon.  Siems  et  al.30  and  Spence31  have  noted  that  d  decreases  as  the 
ribbon  approaches  the  surface  due  to  the  reduced  elastic  repulsion  between  the  partials 
comprising  the  ribbon.  They  have  obtained  an  analytical  expression  for  the  stress  field  of  a 
screw  dislocation  parallel  to  a  free  surface  which  takes  into  account  both  the  surface  effects 
and  anisotropy.  The  stacking  fault  energy,  YSFE.  can  be  corrected  for  surface  proximity 
effects30  using  the  factor,  [1-/1,  given  by  eq’s.  (II.  1-3) 

YSFE  (z  -» oo)  =  YSFE  (z)[l-/l 

where 

[\-f]  =  [1  - (1  +  4z253/x2)-1] 
and 

53  =  (di  -  Ci2)/2C44. 

The  stress  field  at  a  distance  from  the  partial  dislocation,  denoted  by  x,  is  taken  as  the 
observed  separation  between  the  pair  of  partial  dislocations,  d  (=  12-25  nm  in  Fig.  7).  z  is 
the  distance  of  the  dislocation  from  the  free  surface,  and  83  takes  into  account  the 
anisotropy  of  the  material.  In  order  to  correct  YSFE  in  the  STM  experiment  for  surface 
effects,  we  take  z  as  equal  to  either  one  or  two  lattice  interplan ar  c-axis  spacings  (0.335  nm 
or  0.67  nm,  respectively).  These  values  were  chosen  since  the  STM  images  reflect  the 
structure  of  the  first  few  layers  of  graphite.  Using  Eqs.  (II.  1-II.3)  and  the  numerical  result 
from  Appendix  I,  we  obtain  y$fe  =  0. 14  -  2.6  mJ/m2. 

It  is  interesting  to  note  that  if  one  assumes  that  the  value  of  YSFE  obtained  from 
electron  microscopy  studies  (1.5  mJ/m2)  is  the  correct  value  for  the  bulk  stacking  fault 
energy  (i.e.,  as  z  -»  «*>),  the  separation  distance  between  the  partial  dislocations  found  in 
the  STM  experiment  (12-25  nm)  corresponds  to  a  depth  from  the  surface,  z,  of  2.3-7.5  A. 


(n.i) 

(n.2) 

(D.3) 
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This  is  consistent  with  our  assumption  that  the  stacking  fault  network  is  located  within  the 
first  few  layers  of  the  surface. 
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FIGURE  CAPTIONS 

Figure  1.  (a)  Schematic  diagram  showing  a  basal  dislocation  in  graphite  between 
hexagonal  and  i  hombohedral  graphite,  (b)  Projection  of  a  60°  basal  plane  dislocation  on 
the  c-plane.  Two  planes  are  shown:  the  solid  and  dashed  hexagons  are  in  adjacent  planes. 

Figure  2.  STM  image  of  a  (a)  400  nm  x  400  nm  area  of  a  triangular  array  on  the  (0001) 
plane  of  HOPG.  The  array  consists  of  alternating  faulted  (rhombohedral  graphite,  light 
triangular  regions),  and  unfaulted  (hexagonal,  dark  triangles).  Sample  bias  voltage, 
tunneling  cuirent,  and  scan  rate  were  0.24  V,  0.37  nA,  and  26  Hz,  respectively,  (b)  800 
nm  x  800  nm  area  to  the  right  of  the  array  in  (a)  showing  a  row  of  prismatic  dislocation 
loops  which  form  the  boundary  between  the  triangular  array  and  rows  of  ribbons.  Sample 
bias  voltage,  tunneling  current,  and  scan  rate  were  0.12  V,  0.6  nA,  and  26  Hz, 
respectively,  (c)  1200  nm  x  1200  nm  area  to  the  right  of  the  triangular  array  showing  a 
large  area  covered  by  rows  of  ribbons.  The  ribbons  terminate  at  a  defect  site  on  the 
graphic.  Sample  bias  voltage,  tunneling  current,  and  scan  rate  were  0.10  V,  C.78  nA,  and 
26  Hz,  respectively. 

Figure  3.  (a)  STM  image  of  a  5  nm  x  5  nm  area  of  a  dark  triangular  region  (unfaulted) 
showing  the  hexagonal  arrangement  of  the  graphite  lattice.  Bias  voltage  is  0.12  V  and 
tunneling  current  is  0.7  nA.  Inset:  fast  fourier  transform  (FFT)  of  the  image  showing  six 
bright  spots  characteristic  of  hexagonal  graphite.  Periodicities  are  2.6  A.  (b)  STM  image 
of  a  5  nm  x  5  nm  area  near  the  boundary  between  the  lighter  and  darker  triangular  regions. 
The  image  shows  glide  of  the  basal  plane  by  -1  unit  cell.  Bias  voltage  is  0.25  V  and 
tunneling  current  is  1.7  nA.  Inset:  FFT  of  the  image  in  (b)  showing  four  bright  spots  in  a 
rectangular  pattern.  Periodicities  are  2.4  A  and  2.8  A. 
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Figure  4.  STM  images  showing  the  voltage-induced  transition  of  the  star  network  to  a 
triangular  network  as  the  bias  voltage  is  changed  from  0. 10  V  to  0.24  V.  (The  white  dot  in 
the  center  of  Fig.  4(a)  can  be  used  to  identify  the  coordinates  of  the  triangles  in  the 
sequence  of  images.)  (a)  500  x  500  nm  image  of  the  star  network  obtained  at  0. 10  V.  The 
star  network  is  comprised  of  alternating  small  and  large  triangles  pointing  at  180°  from  one 
another,  (b)  300  nm  x  300  nm  image  of  the  same  area  in  (a),  (c)  300  nm  x  300  nm  image 
of  the  same  area  as  in  Fig.  4  (b)  showing  the  reappearance  of  the  triangular  network  -1 
minute  after  changing  the  bias  voltage  from  0.10  V  to  0.24  V.  The  triangles  are  no  longer 
joined  together  at  their  vertices  and  the  base  of  the  triangles  are  trapezoidal  shaped,  (d) 
300  nm  x  300  nm  image  of  the  same  region  taken  1  minute  later.  The  vertices  of  the 
triangles  have  intersected  to  form  bright  points.  Tunneling  currents  are  2.4  nA  and  scan 
rates  are  26  Hz. 

Figure  5.  A  125  nm  x  125  nm  STM  image  of  the  base  region  of  a  triangle  showing  the 
trapezoidal  area  of  light  contrast.  Bias  voltage  is  0.24  V,  tunneling  current  is  2.4  nA  and 
scan  rate  is  26  Hz. 

Figure  6.-  Schematic  showing  the  transition  between  the  triangular-  and  star-shaped 
networks,  (a)  A  triangular-network  of  extended  P  and  contracted  K  nodes,  (b) 
Intermediate  structure  resulting  from  motion  of  one  set  of  parallel  partial  dislocations 
(approximately  left  to  right  across  the  drawing).  The  intermediate  structure  has  a 
trapezoidal  band  at  the  base  of  each  bright  triangle,  (c)  Star-shaped  network  of  fully 
extended  P  and  K  nodes. 

Figure  7.  A  250  x  250  nm  STM  image  showing  pairs  of  partial  dislocations  enclosing  a 
ribbon  of  rhombohedral  (faulted)  graphite.  The  intra-pair  separation  distance  is  marked  by 
the  arrows. 
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